Isoscalar-vector interaction and hybrid quark core in massive neutron stars 
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The hadron-quark phase transition in the core of massive neutron stars is studied with a newly 
constructed two-phase model. For nuclear matter, a nonlinear Walecka type model with general 
nucleon-meson and meson-meson couplings, recently calibrated by Steiner, Hemper and Fischer, is 
taken. For quark matter, a modified Polyakov-Nambu-Jona-Lasinio (mPNJL) model, which gives 
consistent results with lattice QCD data, is used. Most importantly, we introduce an isoscalar-vector 
interaction in the description of quark matter, and we study its influence on the hadron-quark phase 
transition in the interior of massive neutron stars. With the constraints of neutron star observations, 
our calculation shows that the isoscalar-vector interaction between quarks is indispensable if massive 
hybrids star exist in the universe, and its strength determines the onset density of quark matter, 
as well as the mass-radius relations of hybrid stars. Furthermore, as a connection with heavy-ion- 
collision experiments we give some discussions about the strength of isoscalar-vector interaction and 
its effect on the signals of hadron-quark phase transition in heavy-ion collisions, in the energy range 
of the NICA at JINR-Dubna and FAIR at GSI-Darmstadt facilities. 



PACS numbers: 26.60.Kp,21.65.Qr, 97.60.Jd 

I. INTRODUCTION 

The Equation of State (EoS) of neutron star matter is 
closely associated with particles that appear in neutron 
stars. Since the matter in the core of neutron stars is 
possibly compressed to several times of the saturation nu- 
clear density, new particles and even hadron-quark phase 
transitions may appear in the interior of these compact 
objects [1-13]. On the other hand, the EoS of neutron 
star matter is crucial for the macroscopic features of neu- 
tron stars. Each EoS corresponds to one unique mass- 
radius relation of neutron stars by solving the Tolman- 
Oppenheimer-Volkoff (TOV) equation [14]. 

Interestingly, more and more neutron star observa- 
tions, especially, the accurate measurement of the pul- 
sar J1614-2230 with the mass 1.97 ± O.O4M [15], the 
astrophysical observations of X-ray bursts [16-19] and 
thermal emissions from quiescent low-mass X-ray bina- 
ries (LMXBs) in the globular clusters [20-22], gradually 
provide a reliable constraint on the mass-radius relations 
which is tightly connected to the EoS of neutron star 
matter. The analysis of these astrophysical observations 
shows that the radius of a 1.4 solar mass neutron star lies 
between 10.4 and 12.9 km, independent of assumptions 
about the composition of the core [23, 24]. The rela- 
tively small radius of 1.4 solar mass neutron stars means 
the EoS near the suturation density is soft, and the dis- 
covery of massive neutron star J1614-2230 requires the 
EoS is stiff at high densities. The combination of these 
constraints rules out many EoSs of hadron models. Be- 



sides, experimental information from heavy-ion collisions 
(HIC) [25, 26] and lattice QCD simulations [27-33] are 
also available to put some constraints on the EoSs of nu- 
clear and quark matter. 

All these progresses on astrophysical observations and 
laboratory nuclear experiments promote scientists to ex- 
plore the relevant physics behind. The hadron-quark 
phase transition is one of the most concerned topics, and 
so far it is still controversial whether quarks can appear in 
cold neutron star [34-38] . It is also an important topic in 
heavy-ion collisions, and related experiments at medium 
and high densities will be performed in the near future 
on the updated facilities of NICA at JINR-Dubna and 
FAIR at GSI-Darmstadt. 

The hybrid neutron star picture, with direct quark con- 
tributions in the inner core, seems to have problems in 
describing large mass neutron-stars. All that is due to a 
lack of repulsion at high baryon densities of the present 
quark matter effective interactions [39]. This is in fact 
the main point which has motivated the present paper. 

Hadron models based on the relativistic mean- 
field (RMF) theory are usually used to study the prop- 
erties of finite nuclei and nuclear matter. But in litera- 
ture with the same parameter sets, the obtained mass- 
radius (M — R) relations with or without a hadron- 
quark phase transition (e.g., [11, 40-43]) arc not sup- 
ported by the recent analysis of neutron star observa- 
tions [23, 24]. To solve the problem of inconsistency 
between neutron star observations and nuclear matter 
EoS in RMF theory, in this study we will investigate 
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the hadron-quark phase transition with a newly con- 
structed two-phase model. To describe nuclear mat- 
ter, we take the extended Walccka model with general 
nucleon-meson and meson-meson couplings recently cal- 
ibrated by Steiner, Hemper and Fischer [44], which de- 
scribes well the properties of nuclei and nuclear matter, 
and supports recent simulations of the supernovae dy- 
namics. For quark matter, we take the mPNJL quark 
model [41, 42, 45] which shares with QCD global sym- 
metries and the phenomenon of chiral symmetry breaking 
as well as an effective (de) confinement at finite densities 
and temperatures. 

Most importantly, in the description of quark matter, 
we focus on the inclusion of the isoscalar- vector interac- 
tion and its influence on the hadron-quark phase tran- 
sition in the interior of massive neutron stars. With 
the constraints of neutron star observations, our calcula- 
tion shows that the isoscalar- vector channel interaction is 
needed if massive hybrid stars exist in the universe, and 
its strength determines the onset density of quark mat- 
ter, as well as the mass-radius relation of hybrid stars. 
Our previous study [46] also illustrates that this chan- 
nel interaction affects the hadron-quark phase transition 
in heavy-ion collisions at finite temperatures and moder- 
ate densities. Therefore, as a connection with heavy-ion- 
collision experiments, we further give some discussions 
about the strength of isoscalar-vector interaction and its 
influence on the possible phase transition signals from 
asymmetric nuclear matter to quark matter in heavy-ion 
collision experiments. 

The paper is organized as follows. In Section II, we 



describe briefly the two-phase model and give the rele- 
vant formulas of the extended nonlinear Walccka model 
and the mPNJL model. In Section III, we present the 
numerical results, and give some discussions about the 
phase transition in massive neutron stars, as well as the 
connection with the phase transition in heavy-ion colli- 
sions. Finally, a summary is given in Section IV. 



II. THE MODELS 

In the two-phase model, the pure hadronic phase and 
quark phase are described by the nonlinear Walecka type 
model and the mPNJL model, respectively. As for the 
coexisted phase between the pure hadronic phase and 
quark phase, the two phases are connected through the 
Gibbs conditions with the thermal, chemical and mechan- 
ical equilibriums, as well as the global charge neutrality 
condition [47]. 



A. The hadronic model 

Recently one new equation of state of nuclear matter, 
labeled SFHO, was constructed based on the extended 
non-linear Walecka model in RMF theory, and it was 
taken to simulate core-collapse supernova [44]. The ob- 
tained results satisfy the requirements of nuclear physics 
and match well the astrophysical observations. The La- 
grangian of this model is written as 
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where uj pu = d p tu u - d„uj pi p pu = d p p v - d v p^ and 
the scalar meson potential V(a) = (g a a) 3 + \c (g a a) 4 . 
/(er, LO^ujft) takes the general form as 

6 3 

f=Y / w i +Y, b i( uj ^y- ( 2 ) 

i=l j=l 

It was first introduced in [48] to provide additional free- 
dom in varying the symmetry energy, and the parameters 
have been recalibrated recently to fulfill the constraints 
of nuclear physics and astrophysical observations. Due to 



the uncertainties of hyperon-meson couplings, only pro- 
tons and neutrons are considered in this study. Such 
point will be further discussed in Sect. III. Electrons and 
muons, last term of Eq.(l), are included in the calcula- 
tion in order to keep the charge neutrality of neutron star 
matter. 

Under the mean field approximation, the meson field 
equations can be obtained as 
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are the scalar and baryon densities, respectively. 



The thermodynamical potential of the nucleon-meson 
(6) system is 
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where £*(fc) = ^/c 2 + (M, - .gerc) 2 , p* = p t - g^io - 
9p T 3iP f° r nucleons; E*(k) = ^Jk 2 + to 2 , fi* = pi for 
leptons. 

The corresponding energy density and pressure of nu- 
clear matter can be derived as 
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where fi(k) and fi(k) are the fermion and antifcrmion 
distribution functions: 
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The model parameter set labeled SFHO is used in the 
calculation, which can be fixed by fitting the properties of 
symmetric nuclear matter at saturation nuclear density. 
For the details of model parameters, on can refer [44] 



B. The quark model 

For the quark phase, we take the modified three-flavor 
PNJL model with the Lagrangian 



k=0 



(qXkq) 2 + {qilb^kqf 



-K 



det/(g(l + 75)9) + det/(g(l - 75)9) 



-UmA],^[A],T)-Gv{q^qf 



(13) 



where q denotes the quark fields with three flavors, u, d, 
and s, and three colors; too = diag(m M , to^, m s ) in fla- 
vor space; G and K are the four-point and six-point in- 
teracting constants, respectively. The isoscalar-vector in- 
teraction channel is also included. As shown by Eq.(13), 
the p = component of the isoscalar vector interaction 
corresponds to the density operator (qj°q) 2 , therefore it 
is conceivable to expect that the finite-density environ- 
ment brings a significant contribution to this channel. 
Moreover, it is also well known that standard (P)NJL 
models lead to a chiral restoration transition at unphysi- 
cal baryon densities just above the saturation point [49], 
easily reached in heavy ion collisions, without any evi- 
dence of such effect. The presence of a repulsive vector 
field, as considered in the present work, moves the tran- 
sition to more realistic, higher densities. 

This channel interaction reduces the effective quark 
chemical potential, p = p — 2Gvn q . In this study we will 
focus on its influence on the hadron-quark phase transi- 
tion in dense neutron star matter. 

The covariant derivative in the Lagrangian is defined 
as = — iAp. The gluon background field = 
S^Aq is supposed to be homogeneous and static, with 

A = gAQ^-, where 4p is SU(3) color generators. The 
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effective potential W($[A], &[A], T) is expressed in terms 
of the traced Polyakov loop $ = (Tr c L) /Nq and its con- 
jugate $ = (Tr c L')/Nc- The Polyakov loop L is a matrix 
in color space 



L(x) = 'Pcxp 



i I drAi (x, t) 
o 



(14) 



where /3 = l/T is the inverse of temperature and A 4 = 
iA . 

Different effective potentials were adopted in litera- 
tures [41, 50-52]. The modified chemical dependent one 
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was used in [41, 45] which is a simplification of 
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because the difference between $ and $ is smaller at 
finite chemical potential, and 3> = $ at = 0. In the 
calculation we will take the form given in equation. 16, as 
used in [42]. The related parameters, clq = —1.85, a± = 
-1.44 x 1CT 3 , a 2 = -0.08, a 3 = -0.4, are still taken 
from [41], which can reproduce well the data obtained in 
lattice QCD calculation. 

In the mean-field approximation, quarks can be taken 
as free quasiparticles with constituent masses Mi, and the 
dynamical quark masses (gap equations) are obtained as 
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where <pi stands for quark condensate. 

The thermodynamic potential of quark matter in the 
mean-field level can be derived as 
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minimizing the thermodynamical potential 
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All the thermodynamic quantities relevant to the bulk 
properties of quark matter can be obtained from Oq. 
Particularly, the pressure can be derived with P = 
-(fi Q (T,//j - O°(0,0)). From Eq.(18) it is clear that 
the introduction of the isoscalar vector interaction in the 
PNJL model may give important contributions to the 
pressure of quark matter. Indeed the EOS may become 
much stiffcr at high densities for large Gy values. A sim- 
ilar behavior can be obtained in the MIT bag model if 
such term is included. 

As an effective model, the (P)NJL model is not renor- 
malizable, so a cut-off A is implemented in 3-momentum 
space for divergent integrations. The model parameters: 
A = 603.2 MeV, GA 2 = 1.835, KA 5 = 12.36, m u , d = 5.5 
and m s = 140.7 MeV, determined by fitting f nt M n , rriK 
and m v to their experimental values [53] , are used in the 
calculation. 



Concerning the strength of the isoscalar-vector cou- 
pling Gy, there are no explicit constraints at finite den- 
sity. Some efforts to estimate a possible range of val- 
ues for this coupling are briefly described below. For 
the convenience to compare Gy with the strength of the 
isoscalar-scalar interaction G, and for later discussion, we 
define R v = G v /G. 

A very naive estimation, based on the value taken by 
this ratio in the hadronic sector, would give Ry around 
0.6-0.7 [54]. In Rcf.[46] a vector/scalar coupling ratio 
around Ry — 0.2 was obtained from an evaluation of 
only the Fock contributions of the scalar channels. Val- 
ues in the range 0.25 < Ry < 0.5 are derived by a Fierz 
transformation of an effective one-gluon exchange inter- 
action, with Gy depending on the strength of the 17a (1) 
anomaly in the two-flavor model [55, 56]. However, the 
point is that the coupling strength of the direct term 
cannot be fixed, so the total effect of the isoscalar vec- 
tor interaction is still unknown. Other attempts to esti- 
mate Gy are based on the fit of the vector meson spec- 
trum [57]. However, the relation between the vector cou- 
pling in dense quark matter and the meson spectrum in 
vacuum is expected to be strongly modified by in-medium 
effects [52, 58]. 

Because of the uncertainties discussed above, we will 
treat Ry as a free parameter. Aim of the present work 
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is just to get some relevant information from neutron 
star properties. We notice that suggestions to catch in- 
formation on the isoscalar vector interaction may also 
come from heavy ion reactions at next-generation collid- 
ers, such as NIC A and FAIR [46]. A combined study of 
the two aspects appears as a promising tool to get hints 
on the strength of Ry. 

C. The hadron-quark phase coexistence 

The Gibbs criteria are usually implemented for the 
phase equilibrium of a complex system with more than 
one conservation charge. The Gibbs conditions for the 
phase coexistence within a hadron-quark transition in 
compact star are 

T H = T Q , P H = PQ, (20) 

where [i a are usually chosen with fi n and /i e . Under 
the /3 equilibrium without trapped neutrino, the chemical 
potential of other particles including all baryons, quarks, 
and leptons can be derived by 

Ml = 9i^n - qi^e, (21) 

where gi and qi are the baryon number and electric charge 
number of particle species i, respectively. 

The baryon number density and energy density in the 
mixed phase are composed of two parts with the following 
combinations 

P={1-X)Pb+XP% (22) 

and 

s = (1 - X )e H + xe Q , (23) 

where \ is the volume fraction of quark matter. For the 
coexisted phase, the electric neutrality is fulfilled globally 
with 

qtotai = (1 - x) ql P l + * J2 ql P* = °' ( 24 ) 

i—B.l i=q,l 

III. NUMERICAL RESULTS AND 
DISCUSSIONS 

We present in Fig. 1 the EoSs of neutron star matter 
without and with the hadron-quark phase transition for 
different strength of the isoscalar-vector interaction. For 
each value of Ry = Gy/G, the two solid dots with the 
same color indicate the range of the coexisted phase, and 
the cycle marks the largest pressure that can be reached 
in the core of neutron star by solving the TOV equation. 
With increasing vector strength in the quark sector the 
onset of the transition is moving to higher densities since 
the quark pressure is also increasing. 




FIG. 1: (color on line) EoSs of neutron star matter with- 
out and with a hadron-quark phase transition for different 
isoscalar-vector interaction coupling Gv ■ For each value of 
Rv, the two solid dots with the same color indicate the range 
of the mixed phase, and the cycle marks the largest pressure 
that can be reached in the core of neutron star. 



This figure demonstrates that the isoscalar-vector in- 
teraction of quark matter plays an important role in the 
hadron-quark phase transition. In particular the value of 
Rv is crucial for the EoS of neutron star matter at high 
densities. Besides, it shows that the mixed range can 
be reached only inside massive neutron stars. In fact, for 
the case Ry larger than 0.484, the calculation shows that 
quark matter does not appear in the neutron star core. 




FIG. 2: (color on line) The rescaled EoSs as in Fig. 1 with the 
constraints from heavy-ion collisions and neutron star obser- 
vations. The region labeled 68% (95%) gives the 68% (95%) 
confidence level from eight neutron star observation (Ref. [24] 
and refs. therein). HIC marks the constraints inferred from 
heavy- ion collisions [25]. 

To compare with the data from heavy-ion-collision ex- 
periments and neutron star observations, we re-scaled in 
Fig. 2 the EoSs displayed in Fig. 1. It shows all the 
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EoSs can fulfill these constraints at low and moderate 
densities. However, for the case of hybrid star at high 
densities, Ry > 0.1 appears to be favored by recent neu- 
tron star observations of the mass-radius correlation in 
la contours. In the following we will constrain the value 
of Ry using the M — R relations and the accurately mea- 
sured mass of Pulsar J1614-2230. 

As already stated before, in the present study hyper- 
ons are not included in the hadron sector mostly because 
of the uncertainties on the hyperon-meson couplings in 
the nuclear medium. When hypcrons appear in the nu- 
clear matter, just from a degree of freedom counting we 
can expect a sudden softening of the nuclear EoS. Such 
effect can be largely modified from the hyperon interac- 
tions. Indeed this seems to be the indication emerging 
from Heavy Ion data of Fig. 2: the hadronic EoS ap- 
pears to remain rather stiff at densities between 3 and 
5 poj where hyperons would appear in the medium just 
from chemical potential arguments. However we must 
also note that the high density matter formed in HIC 
will certainly differ from the one of neutron-stars, in par- 
ticular for the lack of chemical equilibrium. In any case, 
mass-radius calculations for hybrid neutron stars seem 
to be more sensitive to the model adopted for the quark 
phase than to the hadronic EOS [39, 59]. 

As a further theoretical study, to investigate the prop- 
erties of massive neutron stars with both hyperons and 
quarks, one can fine-tune the meson-hyperon couplings 
and/or introduce the strange meson mediated interac- 
tion between hyperons in the hadron sector, to adjust 
the stiffness of the hadronic EOS. This will be the object 
of a forthcoming analysis. 

In order to see how the isoscalar-vector interaction 
affects the threshold of the hadron-quark phase transi- 
tion, we display in Fig. 3 the relative fractions of dif- 
ferent species as functions of baryon density for Ry = 
0, 0.2, 0.4, respectively. This figure shows that a stronger 
isoscalar-vector interaction postpones the onset density 
of quark matter, and the central density of the corre- 
sponding hybrid star moves to a higher value. However a 
larger Ry also means that the fraction of quark matter is 
smaller in the core of neutron star. Particularly, if Ry is 
large enough, the onset density of quark matter, in fact in 
the mixed phase, will be larger than the central density 
of the neutron star. In this case, no quarks can appear 
in the core of neutron stars. This clearly demonstrates 
the crucial role that the isoscalar-vector interaction plays 
on the hadron-quark phase transition in massive neutron 
stars. 

In Fig. 4 we plot the mass-radius relations of hybrid 
stars with different Ry. The inner (outer) two con- 
tours represent the la and 2a confidence ranges of the 
M - R relations given in Ref. [23] (Ref. [24]), based on 
six (eight) neutron star observations of the X-ray bursts 
[16-19] and thermal emissions from quiescent low-mass 
X-ray binaries (LMXBs) in the globular clusters [20-22] . 
This figure shows that the EoS of neutron star matter 
with the parameter set of SFHO fulfills well the con- 




FIG. 3: (color on line) Relative fractions of different species 
as functions of baryon density with Ry = 0, 0.2, 0.4, respec- 
tively. The blue vertical lines mark the corresponding central 
densities of hybrid stars. 



straints of neutron star observations. This figure also 
gives us an explicit picture of how the strength of the 
isoscalar-vector interaction influences the macroscopical 
properties of massive hybrid stars. The radio timing ob- 
servations of the binary millisecond pulsar J1614-2230 
with a strong general relativistic Shapiro delay signature, 
implies that the pulsar mass is 1.97 ± O.O4M0 [15]. The 
discovery of this massive pulsar rules out many soft EoSs. 
If the isoscalar-vector interaction of quark matter is not 
included, the maximum mass of hybrid stars is 1.88 Mq, 
less than the known maximum neutron star mass. How- 
ever, with the inclusion of this channel interaction and 
taking Ry = 0.055, the obtained maximum mass of a 
hybrid star can reach the lower mass limit of the pulsar 
J1614-2230. Therefore, Ry > 0.055 is required for the 
existence of massive hybrid stars in the universe, and no 
quarks appear for Ry > 0.484 as shown in Fig. 1. 

Finally, we like to present some further discussions 
about the isoscalar-vector interaction in quark matter. 
Apart the general argument about the relevance of this 
term in hadronic matter at low temperature and increas- 
ing baryon density, we have some specific points from 
results of effective non-perturbative QCD models. 

When including this channel interaction in quark 
model, the value of Ry affects the location and emergence 
of critical points of chiral symmetry restoration [56, 60- 
63] . We note that in all the present effective quark mod- 
els with chiral restoration, like NJL or PNJL, show, at 
zero temperature, the chiral transition takes place at un- 
physical low baryon densities, just above the saturation 
value. The inclusion of the isoscalar-vector term, which 
increases pressure and kinetic chemical potential, will 
shift the transition to more realistic density regions. 

Moreover this channel in the effective quark La- 
grangian also influences the onset densities and the ex- 
pected phase-transition signals from charge asymmetric 
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FIG. 4: (color on line) Mass-radius relations of neutron stars. 
The solid curve is the result without quarks and the dash 
curves are the results with a hadron-quark phase transition 
for different Ry- The inner (outer) two contours show the 
la and 2a confidence ranges of the M — R relations given 
in Ref. [23] (Ref. [24]), based on six (eight) neutron star ob- 
servations of the X-ray bursts [16-19] and thermal emissions 
from quiescent low-mass X-ray binaries (LMXBs) in the glob- 
ular clusters [20-22]. 



nuclear matter to quark matter in the two-phase model 
related to experiments in heavy-ion collisions [46]. New 
data about properties of a mixed phase eventually probed 
in high density regions will provide important informa- 
tion on the strength of the isoscalar-vector term in the 
quark interaction. 

In conclusion although this coupling presently cannot 
be determined from experiments and lattice QCD simu- 
lations, there are some good hints about its existence: 

• Compared with the hadron Walecka model, the 
isoscalar-vector interaction of quark matter plays 
a repulsive role similar to the u) meson, impor- 
tant for the properties of nuclear matter in Quan- 
tum Hadron Dynamics model, in particular at finite 
densities and low temperatures. 

• This channel interaction can be derived from higher 



order Fock (exchange) terms or Fierz transforma- 
tions or fitting the vector meson spectrum [46]. 

• The existence of massive hybrid neutron stars as 
shown in this study. 



IV. SUMMARY 

We have studied the hadron-quark phase transition 
in dense neutron star matter with an improved two- 
phase model. The calculations show that massive hy- 
brid stars possibly exist in the universe. In this respect 
the isocalar-vector interaction between quarks is crucial 
for the hadron-quark phase transition. Its strength de- 
termines whether quarks can appear in the interior of 
neutron stars. 

Although the accurate value of Ry is still not known 
by far, neutron star observations can gradually provide 
some constraints on it. In our previous study about 
the hadron-quark phase transition in heavy-ion collisions, 
we have demonstrated that the inclusion of this channel 
interaction postpones the onset density of quark mat- 
ter. The corresponding phase-transition signals, in par- 
ticular on properties of the mixed phase, in the case of 
charge asymmetric nuclear matter to quark matter will 
be strengthened. Then, based on the phase transition 
features of asymmetric strongly interacting matter, we 
also proposed some suggestions to probe the phase tran- 
sition signals in relevant experiments at the FAIR and 
NIC A facilities [46]. The promising conclusion is that in 
a near future the combination of neutron star observa- 
tions and the energy scan of the phase-transition signals 
at FAIR/NICA may provide us some hints on the value 
of Ry, which is helpful for the understanding of quark 
matter interactions and neutron star structure. 
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